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I. INTRODUCTION 



After many successes, standard radiation/matter dominated Big Bang cosmology was 
found inadequate to provide solutions to a number of problems raised when the model 
was studied in more detail in the light of modern data. These problems include the horizon 
problem, the flatness problem, the magnetic monopole problem, etc. Faced with these issues, 
it was clear that a departure from conventional thinking was required and initial assumptions 
needed to be questioned. Specifically, the cosmological constant that was for many years 
in disfavor, was reintroduced and gave a solution to Einstein's equations with exponentially 
growing scale factor, i.e., inflation [1-4]. This immediately solved the problems hsted above, 
since it allowed the universe to be in thermal equilibrium, diluted monopoles, and flattened 
the curvature. Inflation also allowed quantum fluctuations in the early universe to expand 
to super-horizon sizes. Upon re-entry the fluctuations generate the density perturbations 
[5-8] that led to structure [9, 10]. 

As more cosmological data became available [11-13], more detailed inflationary models 
have become necessary to explain it [14]. In the last two decades many inflationary scenarios 
have been analyzed [15-18], but all have one feature in common: homogeneous isotropic 
expansion. However, before the onset of inflation, a typical region of the universe is anything 
but homogeneous and isotropic. It is possible that an asymmetric feature in some region 
could be stretched out by an asymmetric inflation, but remain imprinted through the end 
of inflation, or that a parameter or fleld initially asymmetrically distributed, is diluted to 
a negligible value by the end of inflation, but with an imprint left on the inflated universe. 
Asymmetric features could also be generated by the phase transition that is responsible 
for the inflation. We will investigate these possibihties in models with homogeneous but 
anisotropic expansions. 

In a previous paper [19] we gave exact solutions of Einstein's equations for cases of a 
universe with planar symmetry. These include a universe with cosmological constant plus 
magnetic flelds, cosmic strings or cosmic domain walls aligned uniformly throughout all 
space. In this paper we give exact results that include non-relativistic matter (dust). We 
also give approximate solutions for w 7^ matter. 

The first year WMAP results [20-22] contain interesting large-scale features which war- 
rant further attention [23], [24]. One glaring observational feature is the suppression of 
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power at large angular scales {9 > 60°), which is reflected most distinctly in the reduction 
of the quadrupole C2. This effect was also seen in the COBE results [11, 13]. After the 
COBE experiment, Monte Carlo studies were used to cast doubt on quadrupole suppres- 
sion [25, 26], suggesting the effect could just be statistical. The WMAP analysis [20-24] 
have arrived at similar conclusions. Nevertheless, interesting physical effects are not ruled 
out, especially since the octupole also appears to be somewhat suppressed. Thus it does not 
seem unreasonable to try to model such behavior by altering the cosmological model from 
the standard big bang plus inflation scenario. Intriguingly, the more precise measurements 
of WMAP also showed that the quadrupole C2 and octupole C3 are aligned. In particular, 
the £ = 2 and 3 powers are found to be concentrated in a plane P inclined about 30° to the 
Galactic plane. In a coordinate system in which the equator is in the plane P, the £ = 2 and 
3 powers are primarily in the m — ±i modes. The axis of this system defines a unique ray 
and supports the idea of power in the axial direction being suppressed relative to the power 
in the orthogonal plane. These effects seem to suggest one (longitudinal) direction may 
have expanded differently from the other two (transverse) directions, where the transverse 
directions describe the equatorial plane P mentioned above. Although this effect once again 
could be explained away as statistical [20-24], a realistic physical model that can explain 
some or all anisotropic effects in the WMAP data would be of interest. 

While there are many ways to approach the issue of global anisotropy of the Universe, 
it would be most satisfying to explain global anisotropy by a simple modification of the 
conventional Friedman- Robertson- Walker (FRW) model. To achieve this, one has to consider 
an energy- momentum tensor which is spatially non-spherical or spontaneously becomes non- 
sphcrical at each point in space-time. Such a situation could occur when defects or magnetic 
fields are present. Magnetic fields [27] and cosmic defects [28] can arise in various ways. 
Moreover, it is known that large scale magnetic fields exist in the universe, perhaps up to 
cosmological scales [27, 29] . These considerations motivate us to focus our attention on the 
effect magnetic fields and defects can have on the expansion of the universe. 

As a modest step toward understanding the form, significance and implications of an 
asymmetric universe, we will modify the standard spherically symmetric FRW cosmology 
to a form with only planar symmetry [30]. Our choice of the energy- momentum tensor 
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will result in non-spherical expansion ^ from an initially spherical symmetric configuration: 
an initial co-moving sphere will evolve into a spheroid that can be either prolate or oblate 
depending on the choice of matter content. For the sake of clarity, we first give some general 
properties of cosmologies with planar symmetry. (The universe looks the same from all 
points but the points all have a preferred axis.) Our first example will be a universe filled 
with dust, uniform magnetic fields and cosmological constant. (Some aspects of cosmic 
magnetic fields have been previously studied; see, e. g. Ref. [31].) This is perhaps the most 
easily motivated, exactly solvable case to consider and it will give us a context in which 
to couch the discussion of other examples with planar symmetry and cases where planar 
symmetry is broken. We then describe a number of other exactly solvable planar symmetric 
cases. 

To set the stage, consider an early epoch in the universe at the onset of cosmic inflation, 
where strong magnetic fields have been produced in a phase transition [32-35]. Assuming 
the magnitude of the magnetic field and vacuum energy (A) densities are initially about 
the same, we will find that eventually A dominates. It was estimated [33] that the initial 
magnetic field energy produced in the electroweak phase transition was within an order of 
magnitude of the critical density. Other phase transitions may have even higher initial field 
values [34, 35] , or high densities of cosmic defects. Hence, it is not unphysical to consider 
a universe with magnetic fields and A of comparable magnitudes. If the magnetic fields are 
aligned in domains, then some degree of infiation is sufficient to push all but one domain 
outside the horizon. (Below we also discuss the cases where there is one or only a few 
domains within the horizon.) 

Finally, we should mention that studies of departure from spherical symmetry and/or 
departure from standard inflationary cosmology is an active and controversial area of study. 
A partial list of topics includes polarization of light from astrophysical objects and related 
phenomena [36-40] , the topology of the universe [41-43] , and low I mode suppression [44, 45] . 
The plan of the paper follows. 

In the second section we review the generahzation of a FRW universe to the case of 
planar symmetry. We display the Christoffel symbols, Ricci tensor, and general form of the 

^ By spherical expansion we will mean homogeneous isotropic expansion, while in this paper we will occa- 
sionally use the term non-spherical expansion to describe expansion that has only planar symmetry. 
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energy-momentum tensor with the corresponding Einstein equations and the general form 
of energy-momentum conservation. This section sets up the basic equations to be solved 
and also contains a discussion of thermodynamics in a planar-symmetric universe. We find 
a natural splitting of the elements of T'^'^ into spherically symmetric and anisotropic pieces. 
This procedure provides key insight needed to find exact solutions for the equations of cosmic 
evolution. In the third section we carry out a general analysis of the features resulting from 
planar symmetry, and give various relations and inequalities based on energy conditions, 
many of which involve the eccentricity of the expansion. Various limiting cases are also 
considered. 

Sections 4, 5 and 6 treat a universe filled with cosmological constant, dust and cither uni- 
form magnetic fields, aligned cosmic strings, or aligned cosmic domain walls, respectively. 
We have given each of these exactly solvable cases a separate section so that we can sys- 
tematically compare and contrast them more easily. Graphics and limiting cases are both 
used for this purpose. The magnetic field and aligned cosmic string cases are qualitatively 
similar, and both are substantially different from the case of domain walls. Section 7 con- 
tains our conclusions and a brief discussion of how one would apply our results to density 
perturbations [46]. An appendix has been included to treat the case of matter with generic 
non-zero choice for w, the parameter that describes the equation of state. These results are 
only approximate and so they have been relegated to the appendix to avoid breaking the 
flow of the discussion of exact results presented in the main body of the paper. 

II. UNIVERSE WITH PLANAR SYMMETRY 

To make the simplest directionally anisotropic universe, we modify the FRW spherical 
symmetry of space-time into plane symmetry. (Cylindrical symmetry is, of course, not 
appropriate since it introduces preferred location of the axis of symmetry.) The most general 
form of a plane-symmetric metric (up to a conformal transformation) is [30] 



where a and b are functions of t and z; the xy-plane is the plane of symmetry. We also 
impose translational symmetry along the z-axis; the functions a and b now depend only on 
t. (Examples of plane-symmetric spaces include space uniformly filled with either uniform 
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magnetic fields, static aligned strings, or static stacked walls, where the defects are at rest 
with respect to the cosmic background frame. This situation with defects is artificial or at 
best contrived, but could perhaps arise in brane world physics where walls could be static 
or walls beyond the horizon could be connected by static strings. We will not pursue these 
details here. Of course, any spherically-symmetric contributions (vacuum energy, matter, 
radiation) can be added without altering the planar symmetry.) For the metric (1), the 
nonzero Christoffel symbols are 

= 1^22 — 1 = be , 

Tqi — — fl) Tos ~ ^' 
which results in the following nonzero-components of the Ricci tensor: 

i?% = -(2a + 6 + 2d^ + 6^), 
R\ = R^^ = -{a + 2a^ + ab), 
= -(b + b^ + 2ab). 

To support a symmetry of space-time, the energy-momentum tensor for the matter has 
to have the same symmetry. In the case of planar symmetry this requires 

(T^) = (87rG)-Miag(e,r/,r/,C). (2) 

Here the energy density ^, transverse rj and longitudinal C tension densities are functions 
only of time. The corresponding Einstein equations are 

+ 2ab = ^, (3) 

d + b + d'^ + db + b'^ =r), (4) 
2d + Sa^ = C- (5) 

We also need the equation expressing covariant conservation of the energy-momentum [a 
direct consequence of Eqs. (3)-(5)]: 

e + 2d(e-77) + 6(e-C)=0. (6) 

We consider now the thermodynamics of cosmological models evolving anisotropically. 
Energy density p and pressure p correspond to the spherically-symmetric part of the energy- 
momentum tensor, ^ = A-|-p-|-^, rj — X — p + fj, C^'^^P + C where we have written tildes 
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TABLE I: The components of the energy momentum (2) for various contributions to the matter. In 
the text, tables and figure captions we use the labels A, w, M, S, and W to represent cosmological 
constant, matter with equation of state p = wp, magnetic fields, strings, and walls, respectively. 



Occasionally we label dust with the symbol 0. 

C V C_ 

vacuum energy (A) A A A 

matter (w) p —wp —wp 

magnetic field (M) e — e e 

strings (S) e e 

walls (W) e e 



on the anisotropic parts of the energy- momentum tensor. As in the isotropic case [49], we 
have 

Tdp/dT^p + p, (7) 
where T is the temperature. Up to an additive constant, the entropy in a volume V is 

S^(p + p)V/T. (8) 

Taking V = Vie'^''+'' we find 

S/S^2a + b + p/{p + p). (9) 

For an adiabatic process, the entropy in a comoving volume is conserved; thus the right-hand 
side of Eq. (9) vanishes, 

p+ (2d + 6)(p + p) = 0. (10) 
For matter with the equation of state p — wp, Eq. (10) gives 

p = p.e-(l+-)(2a+6). (11) 

Eq. (10) expresses covariant conservation of the isotropic part of the energy- momentum. 
Since the total energy-momentum is conserved locally [Eq. (6)], the same holds for its 
anisotropic part, 

|+2d(|-ry) + 6(e-C)=0. (12) 
Eq. (12) will be the key to finding our exact solutions. 



III. GENERAL PROPERTIES 



Before considering specific models for the energy-momentum, we first establish several 
general features of an anisotropic universe described by Eqs. (3)-(6). These results will be 
important for conceptual understanding of solutions and asymptotics bounding the effects 
caused by asymmetry and comparing the results for various types of asymmetric components. 

1. We assume that before anisotropic effects became important, the universe had ex- 
panded isotropically. During the initial phase of anisotropic expansion, when anisotropic 
contributions to the energy- momentum are significant, different tension densities in the lon- 
gitudinal and transverse directions cause comoving spheres to evolve into spheroids. At 
a later phase, when all contributions except for the vacuum energy fade away, longitudi- 
nal and transverse expansion rates become equal and the expansion proceeds isotropically. 
Thus, in this universe, each initial sphere develops an eccentricity; whether the resulting 
spheroid is oblate or prolate depends on which tension dominated during the initial phase 
of deformation. 

2. We assume that initially ^ space is isotropic, Oi = bi, and is expanding isotropically, 
di = 6i > 0. Without loss of generality we set ai — O, which is equivalent to a simple rescaling 
of scale factors e" and e**. For expansion in some direction to change into contraction we 

need to cross the point where the expansion rate in this direction becomes zero. Since the 
energy density is positive, from Eq. (3) it follows that there can be no contraction in the 
transverse direction, d > 0. (See Figs. 1, 9, 17.) 

3. If the transverse tension is always smaller (larger) than the longitudinal tension, then 
initially spherical region of space-time has expanded asymptotically to the shape of an oblate 
(prolate) spheroid. Indeed, if r] < (, then Eqs. (4) and (5) give a — h < —lb? + d6 + . 
Consider two cases: (i) when a < b, we arrive at d — 6 > 0, which upon integration leads 
to a contradiction, d — 6 > 0; (ii) when d > b, we find d — b > 2 dt{b'^ — d^), which is 
allowed. Similarly, if 77 > C, then d > b leads to a contradiction, while d < 6 is allowed. In 
the allowed cases, further integration leads to the result: a > b when rj < (, and a < b when 
r] > C- (See Figs. 5, 6, 13, 14, 21, 22 for examples of different cases.) 

^ The subscript "i" refers to the moment of transition from isotropic to anisotropic dynamics. 
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4. On the other hand, longitudinal contraction is possible. To find a moment t — 



when the magnetic field dominates 6* is positive. It follows that only in the magnctic- 
field-dominating case, can space contract in the longitudinal direction; see Figs. 2, 10, 18 
for examples of the different cases. If the initial magnetic field is sufficiently strong, the 
longitudinal size can become smaller then its initial value (curves with 6 < in Fig. 2). 
However, no matter how strong the initial field is, after a sufficiently long period of time 
contraction turns into expansion following the correct asymptotic behavior. 

5. For all known forms of matter the components of the energy-momentum tensor satisfy 
the dominant energy condition [47, 48] ; in our case it says ^ > 0, ^ > 77, and ^ > C (see Table I 
for examples). Evaluating Eq. (6) at the initial time, we find = —ai{3^i — 2'r)i — (^i); from the 
energy conditions it now follows that the energy density does not increase initially, < 0. 
Let us investigate whether the energy can increase at a later time. For this to happen, we 
need to have ^* = for some t = t*; Eqs. (3) and (6) then give (4?7* — 3^^, — C^, ) = ^^,({^, — (*)• 
This equation has a real solution for d* only when the wall contribution dominates; see 
Table I. However, even in the worst case, with only the wall contribution, the energy 
density cannot increase. Indeed, using ^ = 77, we find dl — Since d > by 2, Eq. (3) 
gives 6* = 0. In 4 we found that, unless magnetic field contribution dominates, 6* > and 
thus for walls b cannot become zero if initially hi > 0. This proves that ^ < in all cases. 

6. Let us prove that the transverse expansion rate has its maximum ai t — ti. First, from 
Eqs. (3) and (5) it follows that di = |(Ci~6) ^^'^ so di < because of the energy conditions in 
5. It follows that d initially decreases with time. Suppose, however, that d starts increasing 
and reaches its initial value di for the first time at the moment t^. Equation (5) then gives 

d^: = |(C* — 6) 6 ^ 'C* ^ C* from 5 leads to d* < 0. This is impossible since to reach 
the first point at which d* = d; wc need to have d^, > 0. We thus conclude d < dj, which 
was to be demonstrated. (See Figs. 1, 9, 17.) 

7. Consider now longitudinal expansion. From Eqs. (3)-(5) we have bi — rji — |(^i + Ci)- 
Examining entries in Table I we see that 6i < except for the case when the wall contribution 

^ When we say that some contribution to the energy-momentum dominates we mean by this that the 
contribution is the largest for all times. 
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dominates. Excluding this case, we repeat the argument in 6 and for the point at which 6* = 
6i we find 6* = 77* — + + — of- Using the result in 6, we have 6* <77* — |(^* + C*) < 0, 
which contradicts 6* > that we need to make = hi. We thus conclude that, except when 
walls contribution dominates, the longitudinal expansion rate has its maximum at t = ti. 
(See Figs. 2, 10, 18.) 

8. Next we derive a bound on the spacial volume. Consider two cases: (i) when a < b, 
Eq. (3), ^ < in 5, and d > in 2 lead to d(l + b/d) + b <0 and so 2d + 6 < 0; (ii) when 
a > b, Eq. (3), < in 5, and 6 > in 4 (which is true unless magnetic field contribution 
dominates) lead to d(l + a/b) + ba/b < and so 2d + 6 < again. In both cases integration 
gives the following bound for the volume of space: 

2a + b<{3^{)^t-ti). (13) 

9. The energy condition ^ > C bounds the eccentricity from above. Indeed, from Eqs. (3) 
and (5) follows a — b = {( — ^ — 2d) /2a. Integrating and using d > from 2 and ^ > C from 
5 we arrive at the following bound: 

e""* < di/d. (14) 

From 6 we have di/d > 1 and so Eq. (14) allows a > b. Solving Eq. (5) asymptotically for 
large t, we find d ~ {(/3)^, which turns Eq. (14) into an asymptotic bound 

e"-''<(6/C)^- (15) 

Since e""** is a convenient variable, we define it to be the "pseudo-eccentricity" ^. 

10. When ^ = C (the vacuum energy with any combination of magnetic fields and strings 
aligned in the same direction), similarly to 9 we find 

e'-'-mt^. (16) 

The standard definition of the eccentricity of an elhpse, with semi- major axis of length A = e"", and 
semi-minor axis of length B = e^, is ^^''^^^ = yj g2(a-b) _ \ -y^g g^j-g interested in spheroids that can 
be either prolate or oblate. If a cross section that is tangent to the symmetry axis of the spheroid is an 
ellipse with axes A along the symmetry axis and B normal to that axis, then either one can be larger. 
This distinction is not contained in the definition of eccentricity, so a more appropriate measure for our 
purposes is the ratio ^ = e°'~^ which we will call the pseudo-eccentricity. 
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11. We consider here the case when neither the magnetic field nor the wall contribution 
dominates. Using b > from 4 and b < di from 7, Eqs. (3) and (4) together with the 
condition ^ > from 5 lead to a — b > {d + b) /b, and integration gives the bound 



e"-^ > (b/di) e»/"' 



a/di— 1 



(17) 



The right-hand side of Eq. (17) does not exceed unity, so a < 6 is allowed. Using now 
asymptotic expressions for d and b (which are obtained from Eqs. (3) and (5)), we find 



2(60^ 



exp 



(C/6)^ 



1 



(18) 



Asymptotically ^ K ^^^^ so the bounds (15) and (18) bracket the pseudo-eccentricity 
as follows: 



(A/6)^ exp [(A/6)^ - ll < e"-^ < (6/A)^ 



(19) 



The two bounds in Eq. (19) do not contradict each other since > A. 

12. When ^ = t] (the vacuum energy with walls) we have d — b = (a + b)/b. Using now 
d < — 16 from 8 and integrating we find e'^"'' < (b/di)^. Asymptotically, comoving spheres 
evolve into prolate ellipsoids. 



e"-''<(A/6)i 



(20) 



13. By 6, the transverse expansion rate has its maximum at t — U. Also, for all physical 
matter contributions, the energy density exceeds the vacuum energy (see Table I). Prom 
Eq. (3) we then have d + 2b > X/d, and thus 



a + 2b> A(3/Ci)2(t -ti). 



(21) 



14. Except when contribution of matter with w > dominates, the longitudinal tension 
exceeds the vacuum energy. Equation (5) together with the condition d < dj from 6 then 



gives 



a>X{3ii)~Ht-ti). 



(22) 



15. When A > 0, we have asymptotics ^) C ~ A and d, & ~ (|A)2. In order to find 
asymptotics when A = 0, we first observe that all quantities have power law behavior: 
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^ ~ C^t~^(, etc. Equation (5) gives — 2, Sa — 1. Prom Eqs. (4) and (5) we then have 
either sy > 1, — Sn — 2 or Sj, < 1, = + 1, — 2sy. Examining entries in Table I, we 
conclude that from = 2 it follows that min(se, Sp) = 2; if there is more than one anisotropic 
component, then is the smallest anisotropic exponent. This results in = = 2, > 1. 
(See Table III for examples.) 

16. Let us find how asymptotics for A > in 15 approach their limiting values. Prom 
Eq. (6) we have 

{50'+ (iA)^ (3 - 2n, - n^)5C = 0, (23) 

where — Sr]/S^ and = ^C/^C- When and are constants (for example, when, 
besides the vacuum energy, there is only one other contribution from Table I) or slowly 
varying functions, Eq. (23) gives 6^ oc e~*/*« with characteristic time 

t^ = (|A)-^(3-2n,-nc)-^ (24) 

Por the cases AM, AS and AW, (lA)^^^ equals |, and 1, respectively. Similarly, from 
Eqs. (3)-(5) we obtain 

((5a)'+(iA)^ {3 - 4:n^)Sa - 2n^5b =0, (25) 
(56)-+(iA)^ (4n^-4n^)5d + (3-2n^ + 2nJ56] =0, (26) 

Solving these equations and keeping only the leading terms, we find Sa oc ie~*/*", Sb oc te~^/^i> 
for the case AM with (|A)5t^ = (|A)^t^ = |, and 5a oc 56 oc for the cases AS 

and AW with (|A)^ta equal to 1 and |, respectively, with (iA)^tj, equal to 1 in both cases. 
(See Pigs. 1, 2, 9, 10, 17, 18.) 

Por the reader's convenience we have collected in Table II many of the general results 
proved in 1-16, conditions of their applicability, and examples of matter content when the 
results can be used. 



IV. MAGNETIC FIELDS 



We are now in a position to extend the analysis of a universe with cosmological constant 
plus magnetic fields of Ref. [19] to the more general case that also includes dust. The analysis 
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TABLE II: Summary of general results concerning the system described by Eqs. (3)-(6). Initially, 
space is assumed to be isotropic, ai = 6i = 0, and expanding isotropically, di = di > 0. In 
derivations of some of the results the energy conditions, ^>0,^>r7,^>C, were assumed. The 
symbols in the last column indicate the matter content types for which the corresponding result 
in the second column is applicable. If there are more than two components, then their arbitrary 
combination is allowed. Also, the result can be extended to matter content which is not indicated 
in the forth column provided that it does not exceed the indicated components. 



No. 


Result 


Coudil ions 


A|")i")lical)le (;0 


2 


d > 




A, M, S, W, w 


3a 


a > b 


7 ^ S 


M, S 


3b 


a <b 


V > C 


W 


4 


b>0 


e + 2r?-C>0 


A, S, W, w 


5 


i<o 




A, M, S, W, w 


6 


d < di 




A, M, S, W, w 


7 


6 < di 


^ - 2r/ + C > 


M, S, w 


8 


max (2a + b) 


C + 2r? - C > 


A, S, W, w 


9 


max (a — b) 




A, M, S, W, w 


10 


a — b 


^ = C 


A, M, S 


11 


min (a — b) 


C-2r? + C>0, e + 2r/-C>0 


S, w 


12 


max (a — b) 


^ = V 


A, W 


13 


min (a + 26) 




A, M, S, W, w 


14 


min a 


C>A 


A, M, S, W 



is similar to the case without dust, but adds to it a new variable and requires the use of 
Eq. (10). We will find exact solutions in this case, and in the cases where magnetic field 
is replaced by strings or walls. Lest the reader become complacent with the ease at which 
these solutions have been found, we will show that if we have instead A + M/S/W + w, the 
resulting differential equations are much more complicated and difficult to solve. Numerical 
techniques are still available by which we will solve these equations and present the results 
graphically. 
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In the case of cosmological constant, magnetic fields and 7^ matter, Eqs. (3), (5) and 
(12) with the corresponding entries from Table I give 



+ 2a6 = A + p + e, 
2d + 3d^ = A — wp + e, 
e + 4de = 0. 



(27) 
(28) 
(29) 



From the conservation of the anisotropic part of the energy-momentum, Eq. (29), we find 



a = ^ hi (ei/e). 
Substituting this result into Eq. (28), we arrive at 



(30) 



ee - ^e^ + 2e^(A - + e) = 0. 



(31) 



Eq. (31) does not explicitly involve the independent variable t. To use this fact, we let e 



be the independent variable and introduce a new dependent variable / = ^e^. After this 
change Eq. (31) becomes 



e/'-^/ + 2e^(A-^/;p + e) = 0; 



(32) 



here the prime means differentiation with respect to e. Solving Eq. (27) for b and substituting 
it into Eq. (10) we find (after some algebra) 



efp'-{l + w)p [|/ + e2(A + p + 6)] =0. 



(33) 



The system of coupled differential equations (32) and (33) can be replaced by the following 
system of coupled integral equations: 



/ = |Ae^ + * (pi + 4ei)eT - 8e^ - 2we 



de pe' 



where 



p = p,{e/eOl^'+'"^iJ 



ip — exp 



1 + (1 + w)piei / de (e/ci 



l + w) dee(A + e)/-^ 



(34) 
(35) 

(36) 



We were not able to find the exact solution to the coupled system of Eqs. (32) and (33) 
for arbitrary w. Instead, in the remainder of this section we derive the exact solution for 



14 



w — and defer derivation of an approximate solution for arbitrary w until Appendix A. 
(Another exactly solvable case w = — 1 is not really a separate case since it can be obtained 
from the solution for w = by setting p = and redefining A.) 
For w = 0, Eqs. (32) and (33) give 



/=|Ae^ + |q^(pi + 4ei)e--8e^ 



p = Pi(e/ei)t[l + F(e/ei)]-^ 



A+ (pi + 4ei)(e/ei)3 - 3e 
A + Pi + ei 



where 



F(e/ei) = |(pi/ei)[l + (A + pi)/ei]5 J^^ dxx''^ jA/e; + (4 + Pi/ei)a;i - 3a;} 
Prom Eqs. (11), (30) and (38) we find 



6 = iln 



A+ (pi + 4ei)(e/ei)4 - 3e 



-iln(e/ei)+ln[l + F(e/ei)] 



(37) 
(38) 

(39) 
(40) 



A + Pi + ei 

Finally, time dependence of the above functions p(e), a(e) and 6(e) can be deduced from the 
function e{t), which is given implicitly by 



t-ti 



de 



(41) 



as it follows from f — ^e^ and Eq. (37). 

Prom a physical point of view, the value of e"~^ is the largest for the most anisotropic case; 
for fixed ei, this is achieved for pi = when e""'' ~ (l + ei/A)5 > 1. The case of infinitely 
large pi corresponds to the most isotropic case when e"^'' = 1. It follows that c"^'' > 1 for 
any pi. A careful inspection of the solution embodied in Eqs. (30), (40) and (41) confirms 
this and shows that the space is oblate and its pseudo-eccentricity monotonically increases 
from its initial value (unity) to its asymptotic value. More magnetic field increases the 
anisotropy; more matter reduces anisotropy, but neither can change an oblate ellipsoid into 
a prolate one. 

In the case w = 0, asymptotics from Appendix A simplify as follows: 

e~ei exp -4(A/3)^(i-ti + T)] , (42) 

a~ (A/3)^(t-ti + r), (43) 

b ~ (A/3)5(i - + r) - i In [1 + (pi + ei)/A] + In [1 + F(0)], (44) 

p~Pi[l + (A + pi)/ei]Ml + i^(0)]~'exp [-(3A)5(i - ii + r)] , (45) 
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where 



The corresponding asymptotic for pseudo-eccentricity is 

e"-^~[l + (pi + ei)/A]^[l + F(0)]-\ 



(46) 



(47) 



This form agrees with the general result expressed in Eq. (16). In addition, the lower bound 
for pseudo-eccentricity can be derived: replacing the expression in the braces in Eq. (39) by 
its bound, A/ei + (1 + p\/e\)xi, and integrating, we find 

_ n -1 

Pi + Ci ' 



^a—b 



> 



Pi 



+ 



Pi + ei Pi + ei 



1 + 



A 



> 1. 



(48) 



Hence, the asymmetric expansion is always oblate in this case. 

We finally note that in the case of zero vacuum energy, the exact solution derived above 
simplifies significantly since integrals in both Eqs. (39) and (41) are elementary functions: 

(l + Pi/eO^ I 108-(10 + pi/ei)' 



F(e/ei) = (pi/ei) 



(4 + Pi/eO' 



(1 + PiAi)^ 
108(e/ei)^ - [4 + pi/ei + 6(e/ei)i 



t-ti ^ 2(3ei)-2 



(e/ei)8 ^4 + pi/ei-3(e/ei)4 
4 + pi/ej + 6(e/ei)il [4 + pi/e; - 3(e/ei)3 



(49) 



(e/6i)t(4 + pi/e02 

(10 + Pi/ei)(l + pi/6i)i 
(4 + Pi/e-.r 



(50) 



The resulting form of the solution is then given by Eqs. (30), (38), (40), (49), (50). Using 
this solution we find the following large-time asymptotics: 



e ~ ei 
a ~ 
b ~ 



Mn 



3 

iln 



o 

5(4 + Pi/ei)^(3ei)^t] \ 
|(4 + pi/ei)5(3ei)5i], 

l(p,/eOHA + p,/eO-\3eOHt-ti)], 



p ^ lt-\ 



(51) 
(52) 
(53) 
(54) 
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FIG. 1: Expansion parameter a as a function of t for the case MAw with X = I, pi = 10, ei = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 
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0. 0.2 0.4 0.6 0.8 1. 1.2 t 

FIG. 2: Expansion parameter 6 as a function of t for the case MAw with X = 1, pi = 10, ej = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 

The pseudo-eccentricity becomes 

e'^-'' ~ 1 + 4ei/pi. (55) 

In the following, we will refer to the parameters a and b as planar and axial expansion 
parameters. Note that in order to generate the figures, we have used large values for the 
initial matter density, pi and the magnetic field energy density ei relative to the cosmological 
constant A in order to make the effects of their contributions to the stress energy tensor 
stand out in the graphics. We will do this throughout the paper, but caution the reader 
that we are not implying these are realistic choices of parameters. A realistic choice of 
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FIG. 3: Matter density p as a function of t for the case MAw with A = 1, pi = 10, ei = 200. Curves 
are for w from —1 to 1 with step 0.2 from top to bottom. 




0. 0.05 0.1 0.15 0.2 0.25 t 



FIG. 4: Magnetic field density e as a function of t for the case MAw with X = 1, pi = 10, ei = 200. 
Curves are for w from —1 to 1 with step 0.2 from bottom to top. 

initial conditions would probably be A, pi, and ej all of the same order of magnitude, but 
in this case we would need to plot small differences of parameters instead of plotting them 
directly. Recall that the 6p/p effects found in the cosmic microwave background density 
perturbations are of order 10~^, so observationally one is typically, but not always, looking 
for small effects. 

Figure 1 gives the expansion parameter a as a function of time in a universe filled with 
aligned magnetic fields, cosmo logical constant, and matter. Each curve corresponds to 
matter with a different equation of state parameter w. 
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FIG. 6: Expansion parameters a and b for the case MA with A = 1, = 0. Curves are for 
ei = 0, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000 from top to bottom. 

Each curve in Figure 2 plots the axial expansion parameter b with time in a universe filled 
with aligned magnetic fields and cosmological constant for matter with a variety of choices 
for w. Comparing Figs. 1 and 2 shows a grows faster than b. This implies the expansion is 
oblate, i. e., an initial spherical region expands to an oblate spheroid. 

Figures 3 and 4 show the decay of the matter density and of the magnetic field energy 
respectively, with time using the same initial parameters that were used to generate Figs. 1 
and 2. Figure 5 shows the behavior of the pseudo-eccentricity with time, again for the same 
initial parameters which were used in the previous figures. 
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0. 10. 20. 30. 40. Pi 

FIG. 7: Asymptotic value of the pseudo-eccentricity for the case MAw with A = 1 as a function 
of Pi and ej. Sets of curves are for e"~^ equal to 20,15,10,5 from top to bottom; the abscissa 
corresponds to e""'' = 1. Curves in each set are for w equal to —0.5, —0.25, 0, 0.25, 0.5 from top to 
bottom. 




FIG. 8: Asymptotic value of the pseudo-eccentricity for the case MAw with A = 1 as a function of 
w and Pi for ej = 200. Curves are for e"~^ from 4 to 22 with step 2 from top to bottom. 

In Fig. 6 we have plotted the axial expansion parameter versus the planar expansion 
parameter. Each curve is for a different value of initial magnetic field energy density. Time 
increases along each curve from left to right. Note that a always increases, but for suffi- 
ciently strong initial magnetic fields, after an initial increase, b reaches a maximum, then 
decreases for a time, reaches a minimum, and then increases thereafter. Figure 7 shows the 
asymptotic values of the pseudo-eccentricity as a function of magnetic field energy density 
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and initial matter density for various fixed values of w. As expected, stronger ei leads to 
higher eccentricity, but increasing pi, the spherically symmetric component of T^i, tends to 
dampen the effect. As with previous figures, the parameters in Figs. 6 and 7 were chosen to 
enhance the visualization, not for physical reasons. Finally, Fig. 8 is a contour plot of the 
asymptotic value of the pseudo-eccentricity (the curves are lines of equal asymptotic value 
of e""'') for a range of initial matter densities and equations of state. 



V. STRINGS 

In a somewhat artificial case of cosmological constant plus strings plus matter, the equa- 
tions to solve are 



+ 2ab^ X + p + e, 

2a + 3d^ = A — wp + e, 
e + 2de = 0. 



(56) 
(57) 
(58) 



From the conservation of the anisotropic part of the energy-momentum, Eq. (58), we find 



a = I In (ei/e). 



(59) 



Proceeding with the analysis in a matter similar to Sec. IV, we arrive at the following system 
of equations: 



e/'-i/ + e'(A-«;p + e)=0, 

efp' - (1 + w)p [f / + le2(A + p + e)] = 0, 



or equivalently 



/ = lAe^ + 2e^ + fq ^ (pi - 2ei)e^ - we^ / de pe'^, 



where 



p = A(6/6i)t(l+-)V^ 



ip = exp 



1 + i(l + w)p,e, / de (6/60^+^(1+-)^;/-^ 



n -1 



1(1 + w) / dee(A + e)/-i 



(60) 
(61) 



(62) 
(63) 

(64) 



As in the magnetic field case, we can find the exact solution to the coupled system of 
equations only for non-relativistic matter {w — 0), which we present below. An approximate 
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solution for arbitrary w together with large t asymptotics are given in Appendix B. Setting 
w — 0, Eqs. (62) and (63) reduce to 



/=|Ae2 + 2e3 + |ei^(pi-2ei)e5, 



A + 3e+ (pi - 2ei)(e/ei) 
A + Pi + ei 



(65) 
(66) 



where 



F{e/ei) = f (Pi/ei)[l + (A + p,)/e,]^ J^^ dxx'^ |A/ei + 3x + {p,/e, - 2)xi} 



(67) 



Equations (11), (59) and (66) give 



A + 3e+(pi-2ei)(e/ei)5 ^ 



\ In (e/ei) + In [1 + F{e/e{)\ . 



(68) 



A + Pi + ei 

The time dependences of the above functions p(e), a(e) and 6(e) are found from the function 
e{t) which is given implicitly by 



t-ti = l I de 



Ue^ + e^^ + k; ^(pi-2ei)e 



(69) 



There is a substantial analogy with the magnetic field case: all the statements in the 
paragraph following Eq. (41) in Sec. IV are also correct for the case of strings if one replaces 
the magnetic field density with the string density. 

In the case w — 0, asymptotics from Appendix B simplify as follows: 



e ~ ei exp 



-2(A/3)2(t-ti + r) 
a~ (A/3)i(i-ii + T), 

b ~ (A/3)5(i - + r) - 1 In [1 + (pi + ei)/A] + In [1 + F(0)], 
p - Pi [1 + (A + pi)/ei]^ [1 + F{Q)]-' exp f-(3A)^(t -U + r) 



(70) 
(71) 
(72) 
(73) 



where 



-Xe^ + e^ + U ^(pi-2ei)e 



de|(lAe^)-- 
The corresponding asymptotic for the pseudo-eccentricity is 

[l + (Pi + 6i)/A]Ml + F(0)]-\ 



(74) 



e ~ 



(75) 
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To find the lower bound for this quantity, we replace the expression in the braces in Eq. (67) 
by its bound, A/ei + (1 + Pi/ei)x2, integrate, and find 

-1 



Pi + ei Pi + Ci 



A 



> 1. 



(76) 



This is the same bound we found for the magnetic field case [Eq. (48)], and we see the 
expansion is again of the oblate form. 

Again, for the case A = 0, the above solution is given in terms of elementary functions. 
We find 



F{e/eO = (pi/6i)^(l) 



1 



— arctanh6'(l) + arctanh6'(e/ei) 



where 



^(1) 9{e/e,) 

_i 

t — ti = €; ^ {|(pi/ei — 2) [arctanh6'(l) — arctanh6'(e/ei)] 
-e{l) + {e/eO"^e{e/eO}, 



9{x) = [l + i (pi/ei - 2) a;^ 



The large-time asymptotics are 

a^f{t), 
b^-f{t), 
p~ Pie--^^*\ 

1 1 
where f{t) — te? + |(pi/ei — 2) In {te?). The pseudo-eccentricity becomes 



e ~ e ^ . 



(77) 



(78) 



(79) 



(80) 
(81) 
(82) 
(83) 



(84) 



Figures 9 and 10 plot a and 6 as a function of time for a range of w values. While Fig. 9 
is qualitatively similar to Fig. 1, Fig. 10 shows only a monotonic increase in b, unlike Fig. 
2. (See 2, 4 in Sec. III.) Figures 11 and 12 show the matter density and string density as 
a function of time. These figures are quahtatively similar to Figs. 3 and 4 for the magnetic 
field case, but the curves in the string case are somewhat more compressed. Figure 13 is the 
plot of the pseudo-eccentricity with time for strings. Figure 14 plots a versus b for a variety 
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FIG. 9: Expansion parameter a as a function of t for the case SAw with A = 1, = 10, e; = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 
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FIG. 10: Expansion parameter 6 as a function of t for the case SAw with A = 1, pi = 10, €[ = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 

of values of initial string density. The monotonic increase of b is again apparent, in contrast 
to the results for magnetic fields shown in Fig. 6. 

Figure 15 gives contours of asymptotic values of the pseudo-eccentricity as a function of 
matter density and string density. The results are similar, but somewhat milder than the 
magnetic case, Fig. 7. Finally, Fig. 16. shows a dependence of the asymptotic value of the 
pseudo-eccentricity on the equation of state. The effect is again similar to, but milder than, 
the magnetic field case. 
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FIG. 11: Matter density p as a function of t for the case SAw with X = 1, pi = 10, ei = 200. Curves 
are for w from —1 to 1 with step 0.2 from top to bottom. 




FIG. 12: Magnetic field density e as a function of t for the case SAw with A = 1, pi = 10, = 200. 
Curves are for w from —1 to 1 with step 0.2 from bottom to top. 

VI. WALLS 

If we replace magnetic fields or cosmic strings from the previous two sections by a uniform 
stack of cosmic domain walls, we arrive at another solvable model described by the following 
equations: 

d"^ + 2ab = X + p + e, (85) 
2a + 3d^ = A - wp, (86) 
e + be = 0. (87) 
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FIG. 13: Pseudo-eccentricity e"^"* for the case SAw with A = 1, pi = 10, ei = 200. Curves are for 
w from —1 to 1 with step 0.2 from bottom to top. 




FIG. 14: Expansion parameters a and h for the case SA with A = 1, pi = 0. Curves are for 
ei = 0, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000 from top to bottom. 

Again, we were not able to solve the above equations exactly for arbitrary w\ also, it ap- 
pears to be much harder to arrive at a simple and accurate approximation similar to the 
approximations for the cases of magnetic fields and strings (see Appendices A and B for 
details). We reluctantly restrict ourselves only to the case if = 0. Somewhat surprisingly, 
the analysis will need to be substantially different from the magnetic fields and strings cases. 
To proceed, it is convenient to use the substitution a = | In which transforms the 
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FIG. 15: Asymptotic value of the pseudo-eccentricity for the case SAw with A = 1 as a function 
of Pi and Ei. Sets of curves are for e""'' equal to 20,15,10,5 from top to bottom; the abscissa 
corresponds to e"~^ = 1. Curves in each set are for w equal to —0.5, —0.25, 0, 0.25, 0.5 from top to 
bottom. 
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FIG. 16: Asymptotic value of the pseudo-eccentricity for the case SAw with A = 1 as a function of 
w and Pi for ei = 200. Curves are for e"~^ from 4 to 22 with step 2 from top to bottom. 



Riccati Eq. (86) into an easily solvable linear equation u = jXu. This results in 



a = I m ln7, 

1 — a 



(88) 



where 7 = exp 



{X/3)Ht-ti) 



and 



a 



[l + (pi + 60/A]^-l 
[l + (pi + ei)/A]Ul' 
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(89) 



Energy- momentum conservation, Eq. (87), gives 



6 = In (ci/e). 



(90) 



Using Eqs. (11), (85), (88), and (90), we find 

1 + 

7"'-|-(T\i — (7/ I 

where 



l + (7 /7^-<7^' j-, , 3ei l + cr , Pi _^ 7^ - 1 

7^ + (7 



2A(i-a)^ 



ri^(7) + ^^(l-^ 



2A 



(91) 



F(7) 



dx ■ 



(92) 



/l (X3 + (7)2 ■ 

The function ^(7) can be written in terms of the hypergeometric functions, but the expres 
sion is comphcated and not illuminating, so we will not give it here. 

In the case A = 0, the above exact solution simplifies significantly. Instead of Eq. 
we now have 



a = f ln(l -C + C7t), 
where C — 1 — |[3(pi + ei)]^, and Eq. (91) becomes 



(93) 



[1-C + Ct)^ -1 



(94) 



Returning to the case of A 7^ 0, the large t wall energy density is 

,-(A/3)2(t-ti) 



e ~ ei e 



the transverse scale factor is 



1 + a 2A ^ ^ 2A 1 + a 



(95) 



a~ (A/3)2(i-ii) - I In (1 -cr). 



(96) 



the longitudinal scale factor is 

6- (A/3)^(t-ti) +ln 
and the matter density is 



1 + a ^2A^^^2A 1 + a 



(97) 



p ~ e 



1 — (7)3 3ei , ,4 , , Pi (1 — (7)3 



-1 



(98) 
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FIG. 17: Expansion parameter o as a function of t for the case WAw with A = 1, pi = 10, ei = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 

From Eqs. (96) and (97) the pseudo-eccentricity is 

(99) 

In the case of zero vacuum energy, the asymptotics read: 

e ~ const r^, (100) 
a ~ |lnt + const, (101) 
6 ~ 21nt + const, (102) 

10 

p~ const (103) 

As opposed to the cases of magnetic fields and strings, the pseudo-eccentricity vanishes for 
large t, since e*^" ~ )f:~3. 

The plots for walls shows a number of qualitative differences with the previous cases 
of magnetic fields and strings. The expansion parameters a and b change with time (see 
Figs. 17 and 18) more like the string case, where both grow monotonically; but note that 
now b grows faster than a, which is the reverse of the behavior seen for strings. This means 
the expansion for walls is always prolate, unlike the expansions for strings and magnetic 
fields which are always oblate. 

For walls, the matter density. Fig. 19, falls faster with time than its string and magnetic 
field counterparts due to the fact that the overall expansion, and therefore density dilution. 



e ~ 



1 — a Sei , X 2 , Pi (1 — a) 
+ —^1 - a oo + ^- - 
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FIG. 18: Expansion parameter 6 as a function of t for the case WAw with A = 1, pi = 10, ei = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 
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FIG. 19: Matter density p as a function of t for the case WAw with A = 1, = 10, e; = 200. 
Curves are for w from —1 to 1 with step 0.2 from top to bottom. 

is faster for walls. On the other hand, due to the nature of its contribution to the stress 
energy tensor, the wall energy density. Fig. 20, falls more slowly than the magnetic field and 
string energy densities. 

Figure 21 for the pseudo-eccentricity provides another way of visualizing the pro lateness 
of the wall expansion, since e""'' is always less than one in this case. Figure 22 for a versus 
b represents the degree of prolateness for wall expansion for a variety of initial conditions. 
Figure 23 shows asymptotic pseudo-eccentricity contours for walls as a function of initial 
matter and wall energy densities, and Fig. 24 gives asymptotic pseudo-eccentricity contours 
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FIG. 20: Magnetic field density e as a function of t for the case WAw with A = 1, pi = 10, ei = 200. 
Curves are for w from —1 to 1 with step 0.2 from bottom to top. 




FIG. 21: Pseudo-eccentricity ^ for the case WAw with A = 1, pi = 10, ei = 200. Curves are for 
w from —1 to 1 with step 0.2 from top to bottom. 

for walls as a function of matter density and equation of state parameter w. 

VII. CONCLUSIONS 

Einstein's equations for magnetic fields that extend across the Universe have been consid- 
ered elsewhere. Examples include cylindrically symmetric magnetic geons, exact solutions 
for planar geometry with magnetic fields and dust, and asymptotics [31]. None of these 
studies contain exact solutions with cosmological constant and magnetic fields (AM). We 
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FIG. 22: Expansion parameters a and b for the case WA with \ = 1, pi = 0. Curves are for 
ei = 0, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000 from bottom to top. 
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FIG. 23: Asymptotic value of the pseudo-eccentricity for the case WAw with A = 1 as a function of 
Pi and Ei. Sets of curves are for e"~^ equal to 0.02,0.03,0.04,0.05 from top to bottom; the abscissa 
corresponds to e"~^ = 1. Curves in each set are for w equal to —0.5, —0.25,0,0.25,0.5 from top to 
bottom. 

have not only given exact solutions to the AM case, but also have found exact AM plus 
dust solutions. In addition we have exact solutions when the magnetic fields are replaced by 
uniform arrangements of cosmic strings or cosmic domain walls. Finally we have given ap- 
proximate solutions in all these cases where dust can be replaced by matter with an arbitrary 
value of w in its equation of state. All our solutions have planar symmetry. 

The cosmic microwave background (CMB) and other modern cosmological data are of 
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FIG. 24: Asymptotic value of the pseudo-eccentricity for the case WAu) with A = 1 as a function 
of w and p\ for ei = 200. Curves are for e"~^ from 0.0115 to 0.0120 with step 0.0001 from bottom 
to top. 

such high quaUty that it is now possible to study aspects of the Universe that were previously 
completely out of reach. In order to carry out these investigations, it may be necessary to 
go beyond the homogeneous isotropic big bang/inflationary cosmology and compare the 
data with less symmetric but perhaps more realistic models. In the case of planar symmetry 
studied here, an understanding of the density perturbations and structure formation requires 
perturbing around planar sjTiimetric solutions. Here we have taken a step in that direction 
by considering a planar symmetric universe with eccentric expansion, and have shown exact 
solutions can be obtained even when the eccentricity is large. This will allow a density 
perturbation analysis to be carried out in these cosmologies, which in turn can be compared 
with CMB data and galaxy structure and correlation data [46]. 

It is not just a mathematical exercise to consider planar symmetry. We know that mag- 
netic fields and cosmic defects can be produced in the early universe. In the case of magnetic 
fields, their energy density e at its production epoch can be a substantial fraction of the mat- 
ter density p, and this can cause spherical symmetry to be lost in a cosmology. If the typical 
magnetic domain size D is small compared to H^^ , then the local expansion is eccentric, 
while the average global expansion remains spherical, while ii D > H~^, then the whole 
universe expands eccentrically until D comes within the horizon. Also, ii D < initially, 
a period of inflation can push regions of size D outside the horizon, and we are again in a 
situation of eccentric expansion. 
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The planar symmetric cases of cosmic strings and domain walls are somewhat more 
artificial, since they are assumed to be static and aligned. However, this may not be totally 
unrealistic when considered from the perspective of more fundamental theories. For instance, 
certain AdS/CFT theories derived from string theory have parallel walls, and other theories 
with branes can have strings connecting them. If two parallel walls, both outside the horizon, 
were connected by strings, then the strings would be expected to be parallel on average even 
if they had some dynamics. Based on the above remarks, and with the knowledge of the fact 
that aligned walls and strings both produce planar symmetry, we have given exact solutions 
for these cases as well. 

Even though magnetic fields, string, and wall systems all have planar symmetry, the form 
of their energy momentum tensors differ. For magnetic fields, T^^, is tracclcss, and so this 
case has similarities with a radiation filled universe. For strings and walls, the trace of T^i, 
does not vanish, so there are some similarities with the non-relativistic matter component. 
Strings and walls are under tension, so they also have some similarities with vacuum energy. 
To see all these properties, we have solved the equations of motion exactly for many cases 
of interest. The large-time behaviors of these solutions are summarized in Table III. 

In all these cases, the universe undergoes eccentric expansion and in some instances 
eccentric infiation. Our analysis is completely general, and in order to apply these results, 
more input is necessary, e. g., initial conditions need to be specified, perhaps as derived from a 
model with early universe phase transitions. Time scales need to be fixed, e. g., when did the 
phase transition take place? For instance, for magnetic field production, a phase transition 
not far above the electroweak scale may be effective in producing eccentric effects and at the 
same time remaining compatible with other requirements on the cosmological model, e. g., 
successful baryogeneses. If the magnetic field production scale were too high, then there is 
a danger that all the eccentric effects could be washed out. 

As stated above, with exact planar symmetric solutions at hand, we are now in a position 
to begin density perturbations analysis [46]. To apply the results of this paper, it will be 
necessary to consider how the spectrum of density perturbations are effected by asymmetric 
expansion. Since perturbations get laid down by quantum fluctuations and then asymmet- 
rically expanded in our models, any initial spherical perturbation becomes ellipsoidal. After 
a while, the expansion becomes spherically symmetric again, but as long as perturbations 
remain outside the horizon they stay ellipsoidal. Only after they reenter our horizon will 
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TABLE III: Summary of large-time behavior for various quantities for ten different cases of universe 
content. For each choice of an anisotropic component, magnetic fields (M), strings (S) or walls 
(W), matter with u; = or with < < 1 is included and cosmological constant is either present 
(A) or absent. Only the leading terms in asymptotics are given and i = (A/3)2t. 
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they be able to adjust (they will probably start to oscillate between prolate and oblate with 
frequency that depends on size and overdensity) . So if the perturbations are just entering 
at last scattering they should be elhpsoidal. The smaller they are at last scattering, the 
more they have oscillated and if damped, the closer to spherical they should be. Hence, 
the larger scale perturbations (corresponding to smaller £) will have a better memory of the 
eccentric phase. This would appear to agree with what seems to be hinted at in the WMAP 
observations: more distortion of the low i modes. However, a detailed phenomenological 
analysis needs to be carried out to confirm these facts. 

To summarize, what we need are modes that expanded eccentrically to be entering the 
horizon at the time of last scattering and then to feed this information into a Sachs- Wolfe 
type of calculation. This is a most interesting and challenging calculation, since it requires 
a full reanalysis of the density perturbations in eccentric geometry. In this paper we have 
moved toward this goal. We have carried out exact calculations of the evolution of a variety 
of Universes with asymmetric matter content. In some cases, namely when w is neither 
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zero nor minus one, we have been forced to use approximate methods. We have explored 
the asymptotic behavior of both the exact and approximate cases. Our results provide a 
starting point for the analysis of density perturbations in asymmetric cosmologies. WMAP 
and its successors will be able to either bound or detect effects of asymmetric inflation and 
we have taken the first steps in the theoretical exploration in that direction. 

Finally, we make a few comments about the case where there are multiple magnetic 
domains within the cosmological horizon. (A similar discussion would apply to strings and 
walls.) If the domains are randomly oriented then what one should expect is eccentric 
expansion within each domain, with dependence on the local value of the cosmological 
constant, magnetic field strength, and matter content. Locally there is planar symmetry, 
but globally the Universe would look isotropic if averaged over many domains. One effect of 
the averaging would be an alteration of the power spectrum on scales of order of the domain 
size. This assumes the domains have a preferred size, that is probably on order of the horizon 
size when they were produced, if the associated phase transition was second order, or on 
the size of the correlation length at production, if the associated phase transition was first 
order. This is in contrast to the density perturbations produced in inflation that typically 
have a flat power spectrum. One would also expect to sec polarization effects to survive in 
the CMB in an isotropic average of magnetic domains. A detailed analysis of these effects 
would take dedicated numerical studies. 
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APPENDIX A: APPROXIMATE SOLUTIONS FOR A + M + MATTER WITH 
ARBITRARY w 

We develop a simple approximation by expanding around e = 0. Eliminating p in 
Eqs. (32) and (33), we find 



2wf 



e'f" - 'ief + f / + 2e3 



(1 + w) [ef - i(ll - 3w)f + 2(1 + w)e'{X + e)] . (Al) 
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A power series solution in e to this equation is 



O 



,i(15+3«;) 



(A2) 



The approximate solution for p(e) can then be found from Eqs. (35) and (36). [A much faster 
way to calculate p is to use Eq. (32) directly. The result, p fa pi(e/ei)3(^+"'), is unacceptably 
inaccurate as is clear from both the exact solution (38) for ty = and the asymptotic form 
(A6) below for arbitrary w.] The functions a(e) and 6(e) are given by Eqs. (30) and (11). 
Finally, time dependence of the above functions p(e), a(e) and 6(e) can be deduced from the 
function e{t), which is given implicitly by 



t-t,^l£ de {iAe^ + ier^ [pi(e/ei)^ + ^ e^ - e=^} 



(A3) 



as it follows from / = -^e and Eq. (A2). 

Comparing Eqs. (37) and (A2), we notice that the above approximate solution becomes 
exact for w — 0. In addition, being an expansion in small e, the approximate solution gives 
correct asymptotics for large t. To find the behavior of various quantities for large t, we 
need the corresponding asymptotic of the integral in Eq. (A3). When A > 0, the integral 
diverges for small e, and so we extract this divergent part first; this results in 

t-t,^l{3/X)hn(ei/e)-r, (A4) 

where 



Similarly extracting the divergent part of ip for small e, we find 



Ul+w)(l+w)<l> 



where 



= -|(ei/A) + £ dee(A + e) {(fAe^)-' - f-'] , 



F(0) = (1 + w)p,e, / de (e/eO^+t^^+^V/"'- 
Jo 

Finally, this results in the following asymptotics: 

e~eie-^^, 
a ~ cr, 

6~(7-0 + (l + w)-^ In [1 + F(0)] , 

p ~ + F(0)]-^exp [-(1 + w){3a - 0)] , 
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(A7) 
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(AlO) 
(All) 
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where a{t) — {X/3)^{t — ti + r). For large t, both scale factors grow linearly (as in the 
isotropic case driven by the cosmological constant only). Due to anisotropy introduced by 
the magnetic fields, however, the space has expanded more transversally than longitudinally. 
This difi^erence is characterized by the pseudo-eccentricity whose asymptotic form in this case 
is 



e ~ 



e"^ [1 + F(0)]-'/('+"') . (A13) 



When A = 0, the asymptotics depend on the range of the parameter w; in the most 
interesting case, < w < 1, they are: 

e~3-te~3rt, (A14) 
a -fin (eft), (A15) 

4(l-w) ^_2 

Thus, in the absence of constant negative pressure from the cosmological constant, 
anisotropy causes the space to infinitely expand in the transverse directions and infinitely 

contract in the longitudinal direction. This results in pseudo-eccentricity diverging for large 

1 

t: e«-^ ~ (efi)2"'/(i+"'). In the case of dust {w — 0), the asymptotic value for the pseudo- 
eccentricity is finite, in agreement with Eq. (55). 



APPENDIX B: APPROXIMATE SOLUTIONS FOR A S -h MATTER WITH 
ARBITRARY w 

Eliminating p in Eqs. (60) and (61), we find 
.e'f"-lef' + 7f + e' 



2wf- 



ef'-lf + e^(X + e) 
A power series solution to this equation is 



1 + w) [ef - i(7 - 3w)f + (1 + w)e\X + e)] . (Bl) 
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(B2) 



The approximate solution for p can then be found from Eqs. (63) and (64). [As in the previ- 
ous section, a simple expression p Ki pi(e/ei)t(^+"'\ which follows directly from Eq. (32), is a 
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poor approximation.] Similar to the case of magnetic fields, the approximate solution (B2) 
becomes exact ior w — 0. 

Proceeding similarly to Appendix A, we find the following large-time asymptotics: 

e~eie-2^, (B3) 

a ~ cr, (B4) 

6~ (7-0+ (l + w)-Mn[l + F(0)], (B5) 

p ~ a[1 + F(0)]-^ exp [-(1 + w){3a - 0)] , (B6) 

e«-6^e^[l + F(0)]-'/('+^\ (B7) 

where a{t) = (A/3)5(t — ti + t) and 

T = I d^ I (|Ae^)-^ - [iAe^ + + le7^p,{e/e,y^^ - 2eOei] , (B8) 

= -K^i/A) + I £ dee(A + e) { (fAe^)"' - f-'} , (B9) 

F(0) = i(l + w)p,e, r de {e/e^f+'i^^+^HrK (BIO) 
Jo 

In the case of zero vacuum energy, the asymptotics depend on the range of the parameter 
w] in the most interesting case, < w < 1, they are: 

e~r^ (Bll) 

a ~ In (eft), (B12) 

h^-^\n{e^t), (B13) 

p ^i^-^) t-^ (B14) 

1 + w ^ ' 



As in the magnetic field case, the pseudo-eccentricity diverges for large t: e 
(e.5^)(i+3«')/(i+«'). 
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